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Abstract 

\^ The Largest Empty Circle problem seeks the largest circle centered within the con- 

vex hull of a set P of n points in and devoid of points from P. In this paper, we 
introduce a query version of this well-studied problem. In our query version, we are 
required to preprocess P so that when given a query line Q, we can quickly compute 

^ the largest empty circle centered at some point on Q and within the convex hull of 

I— —I We present solutions for two special cases and the general case; all our queries 

^ run in O(logn) time. We restrict the query line to be horizontal in the first spe- 

^ cial case, which we preprocess in 0(na(n) log n) time and space, where a{n) is 

the slow growing inverse of the Ackermann's function. When the query line is re- 
stricted to pass through a fixed point, the second special case, our preprocessing 
takes 0(na(n)'-^("(")) log n) time and space. We use insights from the two special 
cases to solve the general version of the problem with preprocessing time and space 

O in 0{n^ log n) and O(n^) respectively. 
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1 Introduction 



Facilities that pollute their surroundings are necessary evils. Our cities and in- 
dustrial towns need factories, dump grounds, dams, and nuclear power plants. 
While we cannot eliminate them completely, we would like to locate them far 
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away from human dwellings. The same problem arises on the flip side when 
locating, for instance, a school far away from high crime areas and polluting 
facilities. These scenarios have given rise to a well-studied class of problems 
known as Obnoxious Facility Location. The most basic problem in this class 
is the Largest Empty Circle problem, which takes a set of points P and asks 
for the largest circle with its center inside the convex hull of P and devoid of 
points in P. 

In this paper, we study the placement of an obnoxious facility on a region 
that can be modeled by a line. Consider the obnoxious facility location that 
arises in disaster relief, in which planes on a linear flight path must drop 
personnel in the disaster region. They must, however, be dropped far away 
from points within the region that pose imminent threat. We address this 
flavor of problems by formulating a query version in which we are allowed 
to preprocess the disaster region in a reasonable amount of time, and when 
queried with a flight path, we can quickly provide the best place for dropping 
relief personnel. 

More formally, we study a query version of the Largest Empty Circle problem 
in which we are given set of points P = {pi,P2, ■ ■ ■ ,Pn} strictly inside [0, 1]^. 
When given a query line Q, which we sometime parameterize as Q{t), we are 
to compute the largest empty circle (abbreviated as LEC) with its center on Q 
and within the convex hull of P. The Largest Empty Circle Problem can have 
multiple solutions, but conventionally requires us to report only one solution. 
The query version we study can also lead to multiple solutions, so we keep 
the convention and limit our requirement to one empty circle with the largest 
possible radius. For ease of treatment, we simply call this the largest empty 
circle. 

The Largest Empty Circle problem was first studied in the late 70s and early 
eighties. Toussaint [llj gave an O(nlogn) result and showed that it can be 
extended to the case where the center is constrained to lie within a convex 
polygon. Bose and Wang [3] have shown that O(nlogn) time suffices even 
when the center is to be constrained within a simple polygon. Chew and 
Drysdale [5] studied the problem when the Voronoi diagram and convex hull 
are give as part of the input. They showed that the LEC can be computed in 
0{n) time. The query structure we consider, i.e., requiring the center of the 
circle to be on a query line, is not new either. For instance, Bose et al. j2] 
provide the smallest enclosing circle in 0(?7,log?T.) preprocessing time, 0{n) 
space and O(logn) query time. For more information on obnoxious facility 
location problems, the reader is referred to the survey by Paola [1]. 

In this paper, we present solutions for two special cases and the general case. 
It is well known for the classical Largest Empty Circle problem that the center 
of the LEC will either lie on a Voronoi edge or the convex hull. In Lemma 
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2.1 we show that this is true even when given a query hne, implying that the 
center is at an intersection of the query hne with either a Voronoi edge or the 
convex hulL A query hne can intersect up to 0{n) Voronoi edges implying that 
trivial queries based on this insight alone will not suffice in achieving sub-linear 
query times. The key insight that is novel to this paper is the construction of 
a 3D structure consisting of hyperbolic arcs such that when the query line is 
dropped from z = +00, the point at which it lands on the structure defines the 
center of the LEG. The three cases differ in the way we modify the structure 
so that we can, in O(logn) time for all cases, compute the landing point (and 
hence the center of the LEG). This entails finding the upper envelope of a 
set of functions for which we appeal to the concept of Davenport-Schinzel 
sequences [6l [TJ [10] . 

In the ffist special case, we restrict the query line to be horizontal, i.e., of the 
form y = y^, where y^ G [0, 1] is a constant. We preprocess it in 0{na{n) logn) 
time and space, where a{n) is the slow growing inverse of the Ackermann's 
function. When the query line is restricted to pass through a fixed point, the 
second special case, our preprocessing takes 0(?T,a(?T,)'^^"*^"^Mog?T,) time and 
space. In both cases, we project the 3D structure consisting of hyperbolic arcs 
to a 2D plane so that we can find the landing point. 

For the general case, we assume without loss of generality that the query line 
intersects the x-axis. In addition, we assume that for any fixed query line, 
there are at most 3 LEGs and consider more than three occurrences to be 
degenerate. In Sectio n [5] we describe how this degeneracy can be detected and 



removed. In Lemma 4.1, we show that the x-axis can be divided into O(n^) 
intervals such that each is an instance of the "query line through a point" 
case. This easily leads to an 0(n^a(n)*^*-°^"^'' logn) time algorithm, but we use 
persistent data structures [8] to reduce the preprocessing time to 0{n^\ogn) 
with the data structure taking O(ra^) space. 

The paper is organized as follows. In Section [2| we provide some initial insights 
and a high-level framework for solving all cases. In Section [3} we consider the 
special case in which the query line is guaranteed to be horizontal. The second 
special case in which the query line goes through a fixed point is addressed in 
Section |4j Finally, we provide the solution for the general case, i.e., arbitrary 
query lines, in Section [5| 



2 Characteristics of the Solution 



It has been known that the center of the LEG lies either on a Voronoi edge or 
the convex hull of the points in P yjj. This notion holds true in our situation 
also and is captured by the following lemma. Note that since we require the 
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Fig. 1. The Query Line Q missing all Voronoi edges inside the convex hull. The 
intersections of Q with the convex hull are candidates for the center of the LEG. 

center of the LEG to lie within the convex hull, we only consider portions of 
the Voronoi edges within the convex hull. 

Lemma 2.1 The center of the LEC centered on the query line Q will he the 
intersection point of Q either with a Voronoi edge or the convex hull. 



Proof. Let q E Q he the center of an LEC centered on Q. Let p G P be a 
point closest to q. Consider the Voronoi cell enclosing point p. The query line 
must pass through it (or at least touch it at q) because otherwise, q will be in 
some other Voronoi cell and hence closer to some other point in P. Let Q{t*) 
be the point on Q closest to p. The distance from p to Q{t) will be a convex 
function with its minimum at Q{t*) and strictly increasing on either direction 
of t* . Therefore, the point q will be as far away from Q{t*) as it can go within 
the Voronoi cell of p and within the convex hull of P. This will have to either 
be a Voronoi edge or the convex hull of P. □ 



Lemma 2A_ leads us to the straightforward algorithm of checking all intersec- 
tions of Q with the Voronoi edges and the convex hull edges. The convex hull 
intersections cannot be ignored because in some cases, such as the instance 
shown in Figure [l| the query line misses all the Voronoi edges inside the convex 
hull. Note that the convex hull (and the Voronoi diagram) can be constructed 
during preprocessing in 0(?T,logn) time. 

In the preprocessing step, as a consequence of Lemma |2.1[ we consider the 
Voronoi edges and the convex hull as shown in Figure [2j More precisely, we 
don't consider the Voronoi edges in full, but rather only those segments that 
are encompassed by the convex hull. For the sake of convenience, we use the 
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Fig. 2. The convex hull and the internal segments of the Voronoi edges of a set of 
points. 

generic term, Voronoi edges, to refer to these internal segments. Additionally, 
we break the convex hull into segments delimited by points in P that are on 
the convex hull and the points where the Voronoi edges intersect the convex 
hull. To illustrate, the convex hull in Figure [2] is broken into 19 segments. Note 
that the number of such segments will be 0{n). For the sake of uniformity, we 
treat all the 0{n) convex hull segments and the 0{n) Voronoi edges as a single 
set E. Let n G 0{n) be the cardinality of E. We construct a 3D structure H = 
{hi, /i2, ■ ■ ■ , hn}, where each hj is a hyperbolic arc corresponding to an ej G E. 
The hyperbolic arcs are subsequently transformed in a manner that will allow 
us to query in O(logn) time. Algorithm [l] gives us a high-level framework for 
preprocessing P. While the creation of H (line number |2] of Algorithm [T| is 
common to all versions of the problem that we study, the manner in which 
the structure is transformed to a form that can be queried (line number [s] 
of Algorithm [T]) is quite different for each version and is progressively more 
complicated. We explain the construction of H in this section and defer the 
description of the transformations to later Sections. 

Algorithm 1 Framework for preprocessing P for all cases. 

1: Construct E consisting of the convex hull segments and the internal 
Voronoi edges for P. 

2: Construct the 3D structure H consisting of hyperbolic arcs. This is out- 
lined subsequently in Section |2] and illustrated in Figure [3j 

3: Process H so that the landing point of the query line dropped from z = 
+00 can be computed quickly. 



For each Voronoi edge ej G E, we now describe how and why we construct a 
corresponding hyperbolic arc hj directly above it as shown in Figure [3j Let pi 
and p2 be the two points that induce ej. For convenience, we consider ejit) to 
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Fig. 3. Constructing the hyperbolic arc hj corresponding to a single Voronoi edge 



be the parameterized representation of the point {xt,yt) on ej that is t units 
from the intersection of the Une segment piP2 and the (possibly extended) 
Voronoi edge ej. Intuitively, each point {xt,yt,zt) of the hyperbolic arc hj is 
the elevation of the point {xt,yt) G ej in the +z direction to a height zt that 
equals the euclidean 2D distance from {xt,yt) to either pi or p2- Let d{t) be 
the euclidean distance from pi (or equivalently P2) to ej{t). The height of the 
hyperbolic arc Zt corresponding to the point Cjit) on Cj is equal to the distance 
d{t). Therefore, 



where dist{pi,p2) is the euclidean distance between pi and p2- Hence the arc 
is hyperbolic. 

We note that the above description also holds for the convex hull segments in 
E. Each convex hull segment is closest to at most one point because of the way 
we have segmented the convex hull. This closest point induces the hyperbolic 
arc the same way the two points induce the arc for Voronoi edges. Note that 
this hyperbolic arc will be a straight line for those convex hull segments that 
are incident on a point p in the convex hull of P. This is a natural consequence 
of the way we construct the hyperbolic arcs and does not pose a problem to 
our algorithms. 

While Figure [3] shows the construction of the hyperbolic arc on an infinitely 
long Voronoi edge, our problem is restricted to finding the LEG with center 
within the convex hull, we will only consider the (finite) edges in E. Such finite 
edges will only induce hyperbolic arcs rather than the full hyperbolas. Note 
also that each hyperbolic arc can be constructed in 0(1) time if the points 
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that induce the edg^and the extents of the edge (i.e., the t values between 
which the hyperbohc arc is defined) are given to us. 



One way to interpret this structure in hght of Lemma [2^1] is to drop our query 
fine Q from z = +00 onto the hyperbohc arc structure and report the center 
for the LEG corresponding to the point at which it touches some hyperbohc 
arc and projected straight down onto the z = plane. 



3 Horizontal Query Line 



In this section, we assume that the query line Q will be of the form y = yc, 
where yc is some constant, recall that the points in P lie strictly in [0, 1]^. The 
preprocessing is outlined in Algorithm |2} 

Algorithm 2 Preprocessing P for horizontal query lines. 
1: Construct E consisting of the convex hull segments and the internal 

Voronoi edges for P. 
2: Construct the 3D structure H consisting of hyperbolic circs cis outlined in 
Section |2] 

3: Project each hyperbolic arc in H orthographically onto the a; = plane. 
4: Find the upper envelope of the projected hyperbolic arcs using the algo- 
rithm outlined in [TU] . 

In step number [s] of Algorithm |2| each point (x, y, z) in a hyperbolic arc 
hj will be projected (orthographically) onto the point {0,y,z). The ortho- 
graphic projection preserves the hyperbolic nature of the arcs. Therefore, 
each hyperbolic arc hj becomes a hyperbolic arc h'j on the x = plane. Let 
H' = {h[,h'2,--- , /i'„}. Recall that we are only concerned with the first point 
at which the query line, y = yc, touches the structure when dropped from 
z = +00 onto H. It is easy to see that this corresponds to the upper envelope 
curve in H' at y = yc- Hence, we need the upper envelope of H'. Any two 
projected hyperbolas in H' can intersect at a maximum of two points and are 
partially defined owing to the fact that all hyperbolic arcs in H are constructed 
above the [0, 1]^ region. Therefore, we appeal to the definition of Davenport- 



Schinzel sequences and Theorem 3.3 stated by Sharir and Agarwal [IQl and 
restated here in an abridged manner to capture our requirement. (Please refer 
to the excellent exposition by Sharir and Agarwal [10] for more information 
about Davenport-Schinzel sequences and their application in finding lower and 
upper envelopes.) 

^ Of course, the hyperbohc arcs above the convex huh segment are only induced 
by one point. 
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Definition 3.1 fW^ Let A be an alphabet with n characters and s > be an 
integer constant. A sequence U = ai,a2,--- ,am, where each ai G A, is an 
{n, s) Davenport-Schinzel sequence if it satisfies the following conditions: 

(1) Oj+i for each i < m, and 

(2) there do not exist s + 2 indices ii, Z2, • • • , is+2, where ii < 12 < ■ ■ ■ < is+2, 
such that 

^ii * * * * * * ^ 

and a ^h. 

Definition 3.2 WO^ XAn) = ma.xu \U\, where U is an [n, s) Davenport-Schinzel 
sequence. 

Theorem 3.3 [7^ Given a set H' of n partially defined univariate func- 
tions, the sequence of functions forming the upper envelope is a {n, s + 2) Dav- 
enport Schinzel sequence, where s is the number of points at which two func- 
tions h[ G H' and h'2 G H' can meet. It can be computed in 0{\s+i{n)\ogn) 
time. 

Since our hyperbolas meet at 2 points at most, s = 2. Further, the hyperbolas 
are partially defined. Therefore, from Theorem |3.3| and the upper bounds on 
A functions given in [10], we know that A2+i(n) = A3(n) = 0(na(n)), where 
a{n) is the slow-growing functional inverse of the Ackermann's function. This 
directly leads us to the following corollary. 

Corollary 3.4 The upper- envelope of the set of partially defined hyperbolas 
H' , can be computed in Xzin) \ogn = 0{na{n) logra). 

Proof. In H', all the hyperbolas are oriented upward and hence intersect at 
most at 2 points. Further, since they are partially defined, the proof follows 
from Theorem 13.31 □ 

The upper envelope of H' will be a sequence of maximal intervals of y-values 
such that within each interval we have a single function from H' that domi- 
nates. Hence, when we get a query line of the form y = yc, we merely search 
for the interval that contains yc- This takes O(logra) time. Let Cj E E he the 
edge that induced the hyperbola corresponding to the resulting interval. The 
intersection of Q and cj is the required center of the LEG. 

Theorem 3.5 Given a set P with points in [0, 1]^, we can preprocess P in 
0{na{n)\ogn) time and space such that when we are given a query line Q of 
the form y = y^, where yc G [0, 1], we can report the LEG centered on Q in 
O(logn) time. 
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y-axis 




Fig. 4. Computing the dj{9) function for a single edge e G E. Note that the dj{9) 
function shown on the right roughly corresponds to the edge shown on the left. The 
portion of the curve we need is limited to the angular interval [6j, 9j]. 

4 Query Line Through a Fixed Point 



In this section, we study the special case where the query line Q passes through 
a pre-specified point called the pivot. Without loss of generality, we assume 
that the origin is the pivot. We work out the essential details first assuming 
that P G [0, 1]^; this restricts the pivot to the bottom left corner. Subsequently, 
we will show how this can be generalized to P G [— 1, 1]^. 

In this section, we assume that the set of edges E and the corresponding set 
of hyperbolic arcs H are already constructed. We say that a query line Q 
lands on ej E E ai angle 6 if it makes an angle 6 with the x-axis and an 
LEG is centered at the intersection of Q and ej. In other words, when Q is 
dropped from z = +oo, the hyperbolic arc hj corresponding to ej is the first 
hyperbolic arc it lands on. Drawing from this imagery of the query line landing 
on hyperbolic arcs, we use the phrases "landing on the edge Cj" and "landing 
on the hyperbolic arc hj" interchangeably. In case it lands on more than one 
hyperbolic arc, we simply consider one of them. 

Consider an interval [6*1,6*2] such that for all < 9 < 02, Q lands on a 
hyperbolic arc hj (or equivalently on the corresponding Voronoi edge ej). 
Extending our previous definition, we say that Q lands on hj (or Cj) in the 
angular interval [6^1,^2]- 

To aid in computing the edge on which a query line lands, we first compute a 
function dj{6) for each edge Cj taken individually. We omit the subscript when 
it is clear from the context. Intuitively, it is the height at which a query line Q 
making an angle 6 with the x-axis touches the hyperbolic arc hj constructed 
over Cj. More precisely, if is the maximal range of angle such that a 

query line making an angle 9 E [9^, 9^] intersects edge Cj at a point denoted 
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ej{6), then 

0, Otherwise. 



In the above equation, dist{ej{9),pi) is used to denote the 2D euchdean 
distance between ej{9) and pi.The upper envelope of all these functions, 
D{9) = snp j dj{9'fH provides us the radius of the LEG given a query line 
Q making an angle 9 with the x-axis. With some algebraic and geometric 
manipulations, we get dj{9) to be of the form 



d,{9) 



ai9'^ + a29 + as 



where ai, 02, 03, and 04 are constants dependent on the point in P that induces 
the edge e-,-. 

When given a query line Q passing through the origin and making an angle 
9, we intuitively "drop" the query line from z = +00 and consider the first 
hyperbolic arc it touches. This will be the hyperbolic arc with the largest 
d{9) value. Since the query line can have an arbitrary angle 9, we compute 
the upper envelope D{9) of the hyperbolic arcs in the preprocessing step. In 
the query phase, we can search D{9) to find the appropriate hyperbolic arc 
that Q lands on. Finally, we define the landing sequence of hyperbolic arcs, 
denoted L, to be (/i, h, - ' ' 5 h), to be the sequence of hyperbolic arcs that are 
encountered as a query line is swept from to 7i/2. 

Lemma 4.1 Given a set P of n points, each point lying in [0, 1]^, and a query 
line Q that is restricted to pass through the origin, the landing sequence of 
Q is an (0(n),6) Davenport- Schinzel sequence consisting of partially defined 
functions. 



Proof. The upper envelope curve D{9) that is constructed from the individual 
dj{9) curves of the edges in E determines the landing sequence. (Note that 
\E\ = n & 0{n)). The complexity of D{9) in turn depends on the possible 
number of intersections between any two d{9) curves. If we equate the two 



d{9) functions (given in Equation 4.1) and solve the resulting fourth degree 



equation for 9, we will get up to 4 roots. Therefore, theoretically we can have 
at most four angles at which any two d{9) functions can intersect. Also, we 
know that our hyperbolic arcs (and the useful range of the d{9) functions) are 
of limited size. Therefore, the landing sequence will be an {n, 6) Davenport- 
Schinzel sequence consisting of partially defined functions. □ 

^ We use sup(X) to denote the supremum over a set X of real values, which is 
defined to be the smahest real value that is greater than or equal to every x & X. 



10 



Therefore, from [10], we can construct the upper envelope of the d{9) functions 
for all the hyperbolic arcs (and hence the landing sequence) in 0{X5{n) log n) = 
0{na{n)°'^"^c"'^"^ log n) time, where c > is some constant. Since n G 0{n) and 
c<^(«) grows slower than a{n)"^"''' in the asymptotic sense, we can rewrite the 
running time as 0{na{n)^^°'^^^Hogn) and up to an equal amount of space. 
Although Lemma 4.1 indicates that any two d{6) functions intersect at up to 
4 points, we have been unable to realize this in an example. We believe that 
the sequence will in reality be simpler, but we don't have a proof for it. The 
preprocessing steps are outlined in Algorithm [3j 



Algorithm 3 Preprocessing P for query lines through the origin. 
1: Construct the Convex Hull and the Voronoi diagram for P. 
2: Construct hyperbolic arcs, one for each Voronoi edge and for each convex 
hull segment. 

3: Compute the de function for each hyperbolic arc according to Equation 

sn 

4: Compute the upper envelope D{9) of the set of all dg functions using 
the algorithm outlined in [TU]. Additionally, store the angles at which the 
transitions occur in D{9). 



Given the upper envelope of the d{9) curves and the angles at which the 
transitions occur in the upper envelope, we can, in O(logn) time, find the 
exact hyperbolic arc on which a given query line (passing through the origin 
and making an angle 9 with the x-axis) lands. Substituting the angle 9 in 
the d{9) function for the hyperbolic arc, we can get the radius of the largest 
enclosing circle. The intersection of the query line Q with the edge it lands on 
is the center of the LEC. 



Pivot in Arbitrary Location 



So far, we have worked under the assumption that the pivot point, i.e., the 
point through which the query line must pass, is the origin and the points 
are in [0, 1]^. However, this does not ensure the generality of the solution. In 
particular, what if the pivot point needs to be inside the convex hull of PI We 
can address this by placing the points in P in [—1, 1]^. Now, without loss of 
generahty, the pivot can continue to be the origin. We again need to compute 
the d{9) function for each edge in E, but this time, we sweep a ray (starting at 
the origin) about the origin for the entire 2ti radians. This will not affect the 
asymptotic running time adversely because any two d{9) functions will again 
have at most 4 intersections. A consequence of using rays instead of lines for 
computing the d{9) functions is that when the query line makes an angle 9 
with the X-axis, we have to check the upper envelope at 9 and tt + 6*. 
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5 Arbitrary Query Line 



In this section, we consider tlie general version of tlie problem in which the 
query hue can be arbitrary. Recall that we assume that the points are strictly 
in [0,1]^. Any query line that intersects the convex hull of P must also intersect 
two edges of the [0, 1]^ square. We assume that our query line Q intersects the 
edge on the x-axis; we don't compromise on generality because we can repeat 
this preprocessing for other edges without any asymptotic increase in time or 
space. 

Our approach for this general version builds on the special case studied in 
Section |4] where Q passes through the origin. Like before, we construct a set 
H of hyperbolic arcs in 3D space and seek the point at which Q lands. We 
show that the x-axis (between [0, 1]) can be partitioned into 0{n^) maximal 
segments such that query lines intersecting a given segment induce the same 
landing sequence. Note that the upper-envelopes will vary depending on the 
exact point at which Q intersects the x-axis, but we seek maximal segments 
in which the landing sequence will remain unchanged. 

The algorithm described in Section |4] will not suffice because it requires the 
upper envelope, which lends itself to binary searching. The landing sequence 
is merely a sequence of hyperbolic arcs. The angular range in which each arc 
dominates is not included. So it is not possible to determine whether the query 
line intersects the upper envelope before or after the middle element in any 
portion of the landing sequence. Therefore we start with a detour to show 
how the algorithm described in Section |4] can be modified to work with the 
landing sequence and without the full upper envelope. Subsequently, we show 
that the landing sequences changes incrementally as the point at which the 
query line intersects the x-axis moves along the x-axis. This allows us to store 
the incrementally changing landing sequences in a persistent data structure, 
which can be queried in 0(logr;,) time. 



5.1 Querying with the landing sequence 

In this section, we show that the algorithm described in Section |4] can be 
modified to work even when restricted to using the landing sequence L alone 
and not the upper envelope D{6). For binary search to work, we need to ask 
whether the query line at angle 6 intersects the upper envelope before or after 
the middle element of the sequence and recurse either left or right according to 
the response. As mentioned earlier, we cannot answer this question because we 
don't know the angle at which the d{6) function of each hyperbolic arc in the 
landing sequence starts to dominate over the previous element in the landing 
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sequence. We overcome this limitation by storing a little more information that 
neither increases the running time nor the space (in the asymptotic sense). 

From the definition, in an {n, 6) Davenport-Schinzel sequence such as the 
landing sequence, any two elements a and b can occur in order at most three 
times. Therefore the ■ ■ ■ ,a,b,--- consecutive pair can also occur at most 3 
times. We call them the meetings of a and b. Given a pair a and b and a point on 
the X-axis through which our query line passes, we can find the (at most) three 
meeting angles at which a hands over to 6 in 0(1) time by simply equating the 
two d{6) functions of a and b and solving for 6. Of the four roots and taking 
the extremities of the partial functions, at most three will correspond to a 
handing over to b. Therefore, in the preprocessing phase, we store the meeting 
number (i.e., either first, second or third) along with each meeting without any 
increase in the asymptotic running time or space. This allows us to treat the 
three meetings independently. In Algorithm |3| we add the following line after 
Line number [ij "Construct the landing sequence L = (/i,/2,-"" from 
the upper envelope and additionally store the meeting numbers as a separate 
sequence M = (mi, m2, ■ ■ ■ , where rrij is the meeting number of the 

pair Ijjj+i" 

In the query phase, we perform a binary search on L (in conjunction with 
M) by asking the following question: is the angle 6 that the query line makes 
with the X-axis to the left or right of the rn^th meeting of Ijlj+i, where the 
rn^th meeting of Ijlj+i is the central pair in the portion of the landing sequence 
considered in the current recursion. This can be answered in 0(1) time, thereby 
establishing an O(logn) query time. 

5.2 Partitioning the x-axis 

Let Lc be the landing sequence for query lines passing through (c, 0). Along 
the X-axis, we encounter several maximal intervals of the form [xi,X2], < 
3^1 < 2:2 < 1, during which the landing sequence L^ remains unchanged for all 
c G (xi,X2). As we walk along the x-axis from x = to x = 1, the landing 
sequence changes at a finite number of points that we call events. We now ask 
three questions: 

(1) what are the possible types of changes that an event can induce in the 
landing sequence, 

(2) how many such events exist in [0, 1], and 

(3) how do we compute the events? 

We address these questions in a series of lemmas. In order to prove these lem- 
mas, we consider the following situation to be degenerate: four or more LECs 
occur on a single query line. In other words, the degenerate case happens when 
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a ray lands simultaneously on four or more hyperbolic arcs. Subsequently, we 
will briefly show how to detect and accommodate this degeneracy by slightly 
perturbing the points in P. 



We note that similar situations have been studied previously along with obser- 
vations that are similar to ours. Bern et al. [Ij, for instance, study the changes 
in the topology of a 3D scene when viewed along a straight line flight path. 
Suppose we walk along the x-axis starting from the origin toward x = 1 in 
order to track the changes that we encounter in the landing sequence. Let 
x = c be an event point. For some small e > 0, Lc^e 7^ i-ic+e- Let us first 
assume that the change is an insertion of hyperbolic arcs at a single location 
in the sequence. We will see that either one or two arcs are inserted. Note that 
an insertion when walking in one direction is a deletion when walking in the 
opposite direction, hence our assumption does not affect deletions. Suppose 
the sequence of hyperbolic arcs inserted at x = c is X. We can break up Lc-e 
and Lc+e into subparts consisting of hyperbolic arc sequences (denoted by Si 
and S2), and single hyperbolic arcs (denoted by h, hi, and in two possible 
ways as shown below. 

Case A: Lc-e = Si + h + S2 and Lc+e = Si + h + X + h + S2, or 
Case B: Lc-e = Si + hi + h2 + S2 and Lc+e = Si + hi + J + h2 + S2, 
where "+" is the concatenation operator. In other words, X is either inserted 
in the middle of an existing hyperbolic arc (Case A) or at the meeting point of 
two hyperbolic arcs (Case B). Lemma 5.1 provides us the answer to the first 
question: what types of changes can we encounter? 



Lemma 5.1 If we further divide the cases defined above as shown in the table 
below, Case A* and Case B* cannot occur for non- degenerate input. 



Case name 


Description 


Can occur? 


Case Al 


Case A and X = 1 


Yes (See Figure 


5 


) 


Case A2 


Case A and X = 2 


Yes (See Figure 


6 


) 


Case A * 


Case A and X > 2 


No 


Case Bl 


Case B and X = 1 


Yes (See Figure 


1 


) 


Case B* 


Case B and X > 2 


No 



Proof. Cases Al, A2, and Bl are shown schematically in Figures [5} [6} and 
[7| respectively. Note that each figure highlights a small portion of the the 
hyperbolic arcs viewed orthographically as explained in Section |4j We prove 
the impossibility of occurrence of the two cases separately. 
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0, 



Fig. 5. One hyperbolic arc splits another (Case Al). The hyperbolic arc hi, the 
only element of I, is closer than h to the x-axis on which an orthographic viewer 
is "walking" . Hence, hi moves faster and is visible (and enters the upper envelope) 
after x = c. 




Fig. 6. Two hyperbolic arcs inserted in the middle of another hyperbolic arc (Case 
A2). The hyperbolic arcs hi and /12, which make up the set I, are closer than h to 
the X-axis on which an orthographic viewer is "walking". Hence, they moves faster 
and are visible (and enter the upper envelope) after x = c. 

Case A*: In this case, at our event point x = c, three or more hyperbolic arcs 
that were occluded by the hyperbolic arc h (and hence not present in Lc-J 
appear in Lc. This can happen only if the query line (passing through x = c in 
2D) at the appropriate angle, when dropped from z = +00, will land on the 
three (or more) newly appearing hyperbolic arcs and the occluding hyperbolic 
arc h, thereby leading to four or more LECs, the degenerate case. 
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Fig. 7. A single hyperbolic arc is inserted between two other hyperbolic arcs (Case 
Bl). The hyperbolic arc h, the only element in X, is closer than hi and /12 to the 
X-axis on which an orthographic viewer is "walking". Hence, h moves faster and is 
visible (and enters the upper envelope) after x = c. 

Case B*: If this case were to occur, two hyperbolic arcs h'^ and h'2 must be 
inserted between two other hyperbolic arcs hi and /12 that were present in 
Lc_e. At such an event point x = c, the query line that intersects the x-axis 
at X = c and lands on hi and /i2 will also land on the other two hyperbolic 
arcs h'l and h^-, thereby leading to four LECs. □ 



We now turn our attention to the second and third questions: how many 
event points can we encounter on the x-axis, and how do we compute them? A 



natural consequence of Lemma 5.1 is that each event point is defined by either 



two or three hyperbolic arcs. Lemma 5.2 allows us to compute event points 



when we consider just two or three edges in a Voronoi diagram. Therefore, in 
order to construct all the event points, we have to consider all the hyperbolic 
arc subsets of size two and three. 

Lemma 5.2 All event points induced by a subset of hyperbolic arcs containing 
either two or three arcs can be computed in 0(1) time. 



Proof sketch. For the purpose of this proof sketch, we rely on the ortho- 
graphic viewing of the hyperbolic arcs. Figure |5] illustrates the interaction of 
two hyperbolic arcs to induce a change in the landing sequence. With some 
algebraic and geometric manipulations, we find the point x = c where the 
hyperbola /12 just appears. In similar fashion, we can find the event point at 
which the landing sequence changes for cases A2 and Bl also (see Figures [6] 
and|7|). □ 
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Lemma 5.3 The number of the events on the x-axis from the origin to (0, 1) 
is 0{n^) and the ordered sequence of events can he computed in 0{n^\ogn) 
time. 



Proof. Cases Al, A2, and Bl are all induced by the interaction of either two 
or three hyperbolic arcs. Therefore, we consider all possible pairs and triples 
and compute the set of all the events that each pair or triple can induce 
between [0, 1]. The union of all these sets will give us the set of all events. 
Since there are 0{n^) triples and fewer pairs and given Lemma 5.2, we can 
compute the set of all events in O(n^) time. Therefore, the sorted sequence of 
events can be computed in O(n^logn) time. □ 



In Lemma 5.3, we show that the x-axis can be divided into 0{n^) intervals 
such that each is an instance of the "query line through a point" case. This eas- 
ily leads to a naive 0{n^a{n)^'^°'^^^^ logn) time algorithm in which we store the 
0{rt') difference landing sequences, each of length at most 0(na(?T,)'^'^"*^")^ logn) 
In the query phase, we can find the appropriate landing sequence in O(logn) 
time and subsequently search within the landing sequence using the technique 



described in Section 5.1, again, in O(logn) time. 



We can do better using persistent data structures [8] . We start with a landing 
sequence at the origin. The sequence takes at most 0(na(n)'^(°*^"^^ logn) space 
for the origin. It gets updated at most at each of the Oijn?) event points taken 
in sorted sequence. When we calculate each event, we also store the update 
that takes place; each update requires 0(1) space. In particular, we store the 
following information: 



(1) the X value at which the event occurs, 

(2) the sweep angle at which the event occurs, 

(3) insertion or deletion, 

(4) the sequence of hyperbolic arcs that are inserted or deleted, and 

(5) the meeting numbers to ensure that the landing sequences can be searched 
in O(logn) time. 

For the same reason that the size of the landing sequence is 0(na(n)'^*^"("')^ log n) 
at the origin, the size never gets beyond that at subsequent events. Therefore, 
with the above information, we can use persistent data structures |S] to store 
the information in 0{n^ + na{n)^'^^^'^^^ logn) = O(n^) space. For each update, 
we have to perform an 0(log n) search to find the appropriate location to make 
the update. There can also be a constant number of updates in the meeting 
numbers in the sequence. Each could require a binary search, taking O(logn) 
time. Therefore, the running time for the preprocessing phase is O(n^logn). 
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Algorithm 4 Preprocessing P for arbitrary query lines. 
1: Construct the set E of edges from the Convex Hull and the Voronoi dia- 
gram for P. 

2: Construct the hyperbolic arcs H, each corresponding to an edge in E. 

3: Compute the sorted sequence of event points by considering triples and 
pairs of hyperbolic arcs; additionally store the update information along 
with each event point. 

4: Compute Lq, the landing sequence at the origin and store it in a balanced 
linked binary search tree structure [9]. Embed this linked structure into a 
persistent data structure [8]. 

5: For each event point taken in sequence, update the persistent data struc- 
ture according to the update information. 



Algorithm 5 Querying with arbitrary query lines. 
1: Find the x-intercept, q^, of the query line. 

2: Find the event point c to the left of the x-intercept such that (c, g^;] is 

devoid of event points. 
3: Find the binary search tree version corresponding to event point c and 

perform a binary search for the angle 6 that the query line makes with the 

X-axis. 



In the query phase, we search for the event point that just preceded the x 
intercept of the query line; this takes 0(logr;,) time. We lookup the landing 
sequence induced by the event point, which again takes O(logn) time. This 
will also be the landing sequence at the x intercept. Armed with the landing 
sequence, we binary search in O(logn) time to find the exact hyperbolic arc 
on which the query line lands. This leads us to our final result: 

Theorem 5.4 Given a set P of n points in [0, 1]^, we can preprocess it in 
time O(n^logn) and space 0{n^) such that when presented with an arbitrary 
query line Q, we can report the LEC centered on Q and within the convex hull 
of P in time 0(log?2). 

We stated earlier that we consider the case where more than three LECs occur 
on a query line to be degenerate. This is an essential requirement for Lemma 



5.1 It can be detected when we consider pairs and triples of hyperbolic arcs 
to compute the event points. If more than one event occurs at the same point, 
say x = c on the x-axis, and at the same angle, we can conclude that this 
degeneracy has occurred. We can avoid this degeneracy by slightly perturbing 
the points in P that induce the four or more LECs that cause the degeneracy. 
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